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PREFACE 
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(Code  33A).  The  sponsoring  activity  is  the  Naval  Sea  Systems  Command, 
C.  D.  Smith  (NAVSEA  63R)  Director.  Funding  is  provided  under  Program 
Element  No.  62759N,  D.  E.  Porter,  (NAVSEA  63R3)  Manager.  This  effort 
was  performed  under  Subproject  SF  59-552A. 

The  authors  would  like  to  thank  G.  Swope  (3211)  for  providing  the 
symmetric  arrival  angle  representations,  along  with  D.  Kennedy,  (33A5) 
and  J.  Monti  (3331)  for  their  comments  and  suggestions  during  the 
preparation  of  this  memorandum. 
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Intersensor  cross-correlation  can  be  significantly  affected  by 
acoustic  multipath.  An  analytic  model  has  been  developed  to  simulate 
band-limited  cross-correlation  for  an  environment  in  which  the 
propagation  is  represented  as  four  first-order  bottom  bounce  rays. 

The  "four  path"  model  is  intended  to  be  an  aid  in  understanding  how 
cross-correlation  is  affected  by  bottom  bounce  conditions.  This 
memorandum  describes  the  geometry  and  assumptions  of  the  analytic  four 
path  model,  presents  the  supporting  mathematics,  computes  normalized 
cross-correlations  and  compares  computed  normalized  cross-correlations 
with  correlations  generated  by  the  Generic  Sonar  Model.  The  rms  error 
between  the  two  models  was  found  to  be  less  than  0.1  percent  for  the 
three  test  cases  examined. 
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INTRODUCTION 

An  analytic  model  has  been  developed  to  simulate  cross-correlation 
performance  in  a  bottom  bounce  environment.  This  model  is  based  on 
the  following  assumptions:  1)  band-limited  received  signals,  2)  an 
isovelocity,  lossless  ocean,  3)  the  arrival  at  each  receiver  is  the 
sum  of  four  separate  bounce  paths;  bottom,  bottom-surface, 
surface-bottom  and  surface-bottom-surface,  4)  a  perfectly  reflecting 
bottom,  and  5)  a  perfectly  reflecting  surface  with  180  degree  phase 
shift.  These  assumptions,  while  simplistic,  allow  the  correlation 
function  to  be  expressed  in  closed  form. 

This  memorandum  describes  the  model  assumptions,  geometry,  and 
presents  equations  for  cross-correlation  and  ray  path  travel  time. 
Results  of  a  computer  implementation  are  compared  with  cross- 
correlation  functions  generated  by  the  Generic  Sonar  Model  (ref.  1). 

THE  FOUR  PATH  MODEL 

Figure  lisa  horizontal  projection  of  the  model  geometry,  showing 
the  relative  position  of  source  and  receivers.  Note  the  definition  of 
the  position  of  receiver  1,  and  of  the  bearing  angle,  e.  Ly  is  the 
horizontal  projection  of  the  receiver  separation,  and  is  defined  by 


lh  =  [L2  -  az2]i/2  .  (1) 

where  L  is  the  actual  receiver  separation,  and  az  is  the  depth 
difference  between  receivers. 

With  equation  1,  the  true  range  to  each  receiver  may  be  derived  as 
a  function  of  bearing  angle.  It  can  be  shown  that 


rl  = 

[  r2  + 

<V2>2 

-  Ly • rcos  e  ] 1 1 2 

(2) 

r  2  = 

[  r2  + 

(L„/2)2 

+  Ly • rcos  e  ] 1 >  2 

(3) 

where  r  is  the  horizontal 

range 

to  the 

midpoint  of  the 

recei ver 

basel i ne . 

Figure  2  is  a  vertical  projection  of  the  model  geometry,  showing  the 
bottom  bounce  rays  to  each  receiver.  The  source  is  placed  at  depth 
zQ,  with  the  receivers  at  depths  z.  and  z„,  respectively.  The  time 
delay  between  the  arrival  of  the  signal  at  each  receiver  is  defined  by 


T 


i  j 


(4) 


where  is  the  path  length  from  the  source  to  receiver  1,  for  bounce 

path  i ,  and  is  the  path  length  from  the  source  to  receiver  2,  for 
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bounce  path  j.  When  two  subscripts  are  used  it  is  assumed  that  the 
subscript  i  is  associated  with  receiver  1,  and  subscript  j  with 
receiver  2.  With  this  notation  the  intra-path  del  ay  shown  in  Figure  2 
is  denoted  t..,  being  the  delay  between  the  bottom  bounce  path  at 
receiver  1  ana  the  bottom  bounce  path  at  receiver  2. 

The  total  cross-correlation  may  now  be  written  as  the  arithmetic 
sum  of  the  individual  path  correlations.  For  band-limited  signals, 
the  signal  bandwidth  is  defined  as 


B  =  f H  -  f L  (5) 

where  f„  and  f.  are  the  upper  and  lower  cutoff  frequencies.  Similarly, 
the  center  frequency  is  defined  as 


u> 

o 


2ir 


(6) 


The  correlation  of  paths  i  and  j  is  then  written  as 


x  i  j  (  T )  =  AijCOs[o)o(t-t.j)]sin[TTB(t-T.j)]' 


(7) 


The  amplitude  of  each  correlation  is  the  product  of  the  individual 
path  ampl  i  tildes  ,  or 


Ai  j  =  a  1  i a  2  j 


(8) 


Each  of  the  individual  path  amplitudes  is  a  function  of  the  initial 
signal  amplitude,  Ao ,  and  spherical  spreading  loss  only.  The 
amplitude,  in  decibels,  of  the  signal  at  receiver  k  from  bounce  path 
i  ,  is: 


aki <dB ) 


20Log  A  -  20Log 

10  0  10  1 


(9) 


where  l ^  is  the  length  of  the  ith  bounce  path. 

The  path  lengths  are  easily  derived  for  an  isovelocity  ocean  and 
are  functions  of  depth  and  the  tangent  of  the  arrival  angles.  A 
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vertical  projection  of  the  bottom  bounce  path  is  shown  in  Figure  1. 
It  can  be  shown  that  the  tangent  of  the  arrival  angle  at  receiver  i , 
i  s 


tanc|>B.  = 


2d  -  <z„  +  i,) 


(10) 


and  the  total  bottom  bounce  path  length  is 


*Bi  =  [(d-zo>2  +  C-^2 - )2]W2  +  [(d"zi)2  +  - )2]W2  (U) 

tan*B  .  t a n <(> B  . 


where 

d  =  water  column  depth  in  meters 

z  =  source  depth 

z°  =  depth  of  receiver  i 

r?  =  horizontal  range  to  receiver  i 

Similarly,  for  the  bottom-surface  bounce  path, 


tan^Bsi  = 


(12) 


and 


*BSi  =  [  ( d“zo ) 2  +  (-  ■d~Z.°-.)2]1/2  +  [d2  +  (— - - )2]1/2 


t  a  n  <() 


BSi 


tan^Bsi 


+  U-,2  +  ( - - - )2] 


1  /  2 


tan^Bsi 


(13) 


The  surface-bottom  bounce  path, 
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2d  +  (zo  "  Zi> 
tan+SB1  =  - 


(14) 


and 

*SBi  =  t(^02  *  ( - ~ - )2]1/2  +  [d2  +  ( - - - )*]1/2 

tan^sgi  tan(|>SB. 


+  [(d-Zi)2  +  (— -~-Z.i...)2]^2 

t  a  n  4>  $  b  i 


(15) 


The  surface-bottom-surface  path, 


2d  +  (z0  +  zi  ) 

tan*SBSi  ’  -  (16) 

ri 

and 


^SBSi  =  tzo2  +  ( - ~ - )211/2  +  2*td2  +  ( - - - )211/2 

tan^Bsi  ta"*SBSi 


+  [z-2  +  ( - - )2]^2  (17) 

tan^sBsi 

COMPUTER  IMPLEMENTATION 

The  equations  presented  in  the  previous  section  have  been 
implemented  on  the  Code  33  VAX  system  under  the  name 
DRBO : [HAUCK]4PATH .  Each  bounce  path  length,  and  arrival  angle  is 
calculated  by  a  subroutine  named  DEL**,  where  **  indicates  the  type  of 
path.  For  example,  DELB  calculates  the  path  length  for  the  bottom 
bounce  path.  Attenuation  for  each  of  the  eight  separate  travel  paths 
is  then  calculated  by  the  subroutine  ATTN.  Sixteen  inter-path 
correlations  are  calculated  and  summed,  with  the  resulting  time  series 
normalized  by  the  sum  of  the  four  intra-path  energies.  Each 
intra-path  energy,  where  i=j,  is  defined  by 
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N  1 

Ji  [.1k»  ♦ 


(18) 


where  N  is  the  total  number  of  samples  in  the  observation  window.  The 
total  energy  is  then  the  sum  of  the  four  intra-paths,  or 

en  <19> 


The  final  result  is  then  plotted  using  a  general-purpose  plotting 
routine,  ALLPLOT,  and  written  in  a  Generic  Sonar  Model  compatible 
output  file.  Alternatively,  intermediate  values  of  delay  time,  range 
and  phase  may  be  written  to  the  user's  default  output  device. 

Testing  of  the  model  was  accomplished  by  comparison  with  results 
produced  by  the  GSM.  The  GSM  computes  the  cross-correlation  via  an 
FFT  of  the  synthesized  received  cross-spectrum  (ref.  1).  Three  test 
cases  were  examined.  The  parameters  common  to  all  three  cases  are 
defined  in  Table  I. 


Table  I  -  Parameters  Common  to  All  the  Three  Test  Cases 


Sound  Speed  =  1500  m/s 
Water  Depth  =  1000  m 
Source  Depth  =  25  m 

Range  =  1000  m 

(Source  to  center  of  receiver  baseline) 
fL  =  200  Hz 
fH  =  700  Hz 

Receiver  1  depth  =  50  m 


In  the  first  case,  the  source  is  broadside  (0  deg)  to  the  receiver 
baseline,  and  receiver  2  is  at  the  same  depth  as  receiver  1.  (The 
sample  runstream  required  to  run  the  GSM  for  this  case  is  given  in 
Appendix  I.  The  interested  reader  can  observe  in  the  runstream  the 
appropriate  statements  requried  to  model  a  four  path  isovelocity,  and 
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zero  attenuation  medium).  Figure  3  presents  for  this  case  the 
cross-correlations  obtained  by  both  the  analytic  four  path  model  and 
by  the  GSM.  The  results  are  indistinguishable  by  visual  comparison. 

In  the  second  case  the  source  is  again  at  0°  bearing,  while  receiver  2 
is  now  at  a  depth  of  60  m.  Figure  4  shows  the  comparison  of  the  two 
models  for  the  second  case.  In  the  third  case,  the  source  bearing  is 
now  45°.  Figure  5  shows  the  comparison  of  the  two  models  for  the 
thi rd  test  case. 

The  c ros s -cor rel at i on  time  series  produced  by  each  model  were 
compared  numerically  using  the  program  4ERR.  For  a  given  case,  this 
program  computes  the  rms  amplitude  error  between  the  analytic  and  GSM 
correlations.  In  each  of  the  cases  examined,  the  rms  comparison  error 
was  less  than  0.001,  as  might  be  expected  by  visual  inspection  of 
Figures  3-5.  Also  note  that  the  placement  of  the  correlation  peaks, 
in  time,  agrees  to  within  the  resolution  given  by  the  sampling  rates 
used. 


SUMMARY  AND  CONCLUSION 

To  summarize,  this  memorandum  describes  an  analytic  model  designed 
to  predict  the  cros s -cor rel at i on  of  two  receivers  in  a  four  path, 
bottom  bounce  environment.  By  comparing  the  output  of  the  model  with 
the  GSM,  an  FFT-based  computation  and  the  analytic  solution  may  be 
compared.  For  each  of  the  three  test  cases  examined,  no  significant 
differences  between  the  correlations  obtained  by'the  two  models  were 
observed . 

In  conclusion,  comments  addressing  the  limits  of  the  analytic 
model  are  in  order.  The  entire  model  derivation  is  based  on  several 
assumptions,  restricting  the  type  of  propagation,  sound  speed  profile, 
boundary  interaction,  attenuation  and  the  form  of  the  received 
signals.  The  sound  speed  profile  is  assumed  constant.  While  this 
assumption  implies  a  direct  path  between  source  and  receivers,  the 
model  assumes  bottom  bounce  propagation.  The  only  allowed  paths  are 
the  bottom  bounce,  bottom-surface,  surface-bottom,  and  the 
surface-bottom-surface.  The  ocean  bottom  is  modeled  as  perfectly 
rigid,  with  reflection  coefficient  of  unity.  The  ocean  surface  is 
assumed  perfect  pressure  release,  with  reflection  coefficient  of  -1. 
The  ocean  is  modeled  with  no  attenuation,  the  only  losses  caused  by 
spherical  spreading  of  the  transmitted  signal.  Lastly,  the  received 
signal  is  assumed  band-limited,  with  constant  amplitude  across  the 
signalbandwidth. 

By  employing  these  limitations,  it  is  possible  to  derive  an  exact 
solution  for  the  cross-correlation  function.  The  intent  is  not  to 
replace  more  realistic  computer  models,  such  as  the  GSM,  but  to 
provide  a  model  which  may  quickly  predict,  from  a  closed  form  analytic 
representation,  correlation  performance  in  an  idealized,  bottom  bounce 
ocean.  The  model  thus  can  be  utilized  to  provide  physical  insight 
into  the  effects  of  bottom  bounce  multipath  on  the  correlation 
process  . 
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APPENDIX  I 


COMMENT  TABLE 

RUNSTREAM  FOR  4PATH  COMPARISON 
LOSSLESS,  ISOVELOCITY  OCEAN 
EQUAL  DEPTH  RECEIVERS 
EOF 

RANGE  UNITS  =  M 
TIME  UNITS  =  MS 
TIME  MINIMUM  =  -150.0  MS 
BOTTOM  DEPTH  =  1000  M 
OCEAN  SOUND  SPEED  TABLE 


M  M/S 

0.  1500.0 

250.  1500.0 

1000.  1500.0 

EOF 

RADIUS  OF  CURVATURE  =  999999999  KM 
OCEAN  SOUND  SPEED  MODEL  =  LINEAR 
DEPTH  UNITS  =  M 
VELOCITY  UNITS  =  M/S 
SOURCE  DEPTH  =  25  M 
VOLUME  ATTENUATION  MODEL  =  FRQTBL 
VOLUME  ATTENUATION  TABLE  =  0  DB/KM 
SURFACE  REFLECTION  COEFFICIENT  MODEL 
SURFACE  REFLECTION  COEFFICIENT  TABLE 
BOTTOM  REFLECTION  COEFFICIENT  MODEL  = 
BOTTOM  REFLECTION  COEFFICIENT  TABLE  = 


FREQUENCY  MINIMUM  =  200  HZ 
FREQUENCY  MAXIMUM  =  700  HZ 
FREQUENCY  INCREMENT  =  250  HZ 
RANGE  MINIMUM  =  1000.  M 
RANGE  MAXIMUM  =  1000.  M 
RANGE  INCREMENT  =  1000.  M 
SOURCE  LEVEL  TABLE 
HZ  DB 

1.  -300. 

199.  -300. 

200.  100. 

700.  100. 

701.  -300. 

3100.  -300. 

EOF 

FILTER  EQUALIZER  MODEL  =  TABLE 
FILTER  EQUALIZER  TABLE 
HZ  DB 

100.  -300. 

199.  -300. 

200.  0. 

450.  0. 

700.  0. 

701.  -300. 

3100.  -300. 

EOF 


=  TABLE 
=  0  DB 
•  TABLE 
0  DB 
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TRANSMITTER  BEAM  PATTERN  MODEL  =  TABLE 
TRANSMITTER  BEAM  PATTERN  TABLE 


DEG 

DB 

-90. 

0. 

-11. 

0. 

-5. 

-300 

5. 

-300 

11. 

0. 

90. 

0. 

EOF 

TARGET  DEPTH  =  50  M 
E IGENRAY  FILE  =  TE1 
COMPUTE  EIGENRAYS 
SORT  EIGENRAYS 
PRINT  EIGENRAYS 
E IGENRAY  FILE  =  TE2 
TARGET  DEPTH  =  50  M 
COMPUTE  EIGENRAYS 
SORT  EIGENRAYS 
PRINT  EIGENRAYS 

RECEIVER  BEAM  PATTERN  MODEL  =  TABLE 
RECEIVER  BEAM  PATTERN  TABLE 


DEG 

DB 

-90. 

-300. 

-65. 

-300. 

-64.5 

0. 

-63.5 

0. 

-63. 

-300. 

0. 

-300. 

61.5 

-300. 

62.2 

0. 

63.5 

0. 

64. 

-300. 

90. 

-300. 

EOF 

BEARING 

ANGLE  =  0  DEG 

ARRAY  SEPARATION  TABLE  =  100  M 

NOISE  CORRELATION  DISTANCE  =  5  M 

AMBIENT 

NOISE  SPECTRA  MODEL  =  TABLE 

AMBIENT 

NOISE  SPECTRA  TABLE  =  -200  DB 

FREQUENCY  MAXIMUM  =  3100  HZ 
CORRELATION  FILE  =  CT1 

E IGENRAY 

FILE  =  TE1 

RESET  CROSS-CORRELATION  COEFFICIENT 

COMPUTE 

CROSS-CORRELATION  COEFFICIENT 

RECEIVER 

BEAM  PATTERN  TABLE 

DEG 

DB 

-90. 

-300. 

-65. 

-300. 

-64.5 

0. 

-63.5 

0. 

-63. 

-300. 

0. 

-300  . 

61.5 

-300. 
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62.2  0. 

63.5  0. 

64.  -300. 

90.  -300. 

EOF 

EIGENRAY  FILE  =  TE2 

COMPUTE  CROSS-CORRELATION  COEFFICIENT 
END 
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Figure  1  -  Horizontal  Projection  of  Model  Geometry;  bearing  angle 

6  is  used  by  4  PATH  for  range  calculation,  <j>  is  the  angle 
input  to  the  program. 
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Reverse  Blank 


Figure  2  -  Vertical  Projection  of  Model  Geometry,  showing  Bottom 
Bounce  Angles  and  <j>b2. 
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NORMALIZED  CROSS-CORRELATION 
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Figure  3  -  Generic  Sonar  Model  and  Analytic  Four  Path  Model  Cross-correlation 
Comparison  where  Source  Bearing  Angle  is  0°  and  Receivers  are 
at  Equal  Depth. 
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NORMALIZED  CROSS-CORRELATION 
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™  to  Figure  4  -  Generic  Sonar  Model  and  Analytic  Four  Path  Model  Cross-correlation 

n>  ^  Comparison  where  the  Source  Bearing  Angle  is  0° ,  and  Receiver  2  is 

co  co  at  a  Depth  of  60  meters. 
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NORMALIZED  CROSS-CORRELATION 
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Figure  5  -  Generic  Sonar  Model  and  Analytic  Four  Path  Model  Cross-correlation 
Comparison  where  the  Source  Bearing  Angle  is  45°,  and  Receivers  are 
at  Equal  Depth. 
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